Probing the dynamics of dark energy with novel parametrizations 



O 

< 



o 

U 

Oh; 

I 

o 



(N 

> 
o 



X 



Jing-Zhe Ma' and Xin Zhang ''^ 

'Department of Physics, College of Sciences, Northeastern University, Shenyang 110004, China 
^Center for High Energy Physics, Peking University, Beijing 100080, China 

We point out that the CPL parametrization has a problem that the equation of state w(z) diverges in the far 
future, so that this model can only properly describe the past evolution but cannot depict the future evolution. To 
overcome such a difficulty, in this Letter we propose two novel parametrizations for dark energy, the logarithm 
form w(z) = Wo + wi(i2|2i:^ -In 2) and the oscillating form vv(z) = wo + wiC^^i^j^i^ -sin(l)), successfully avoiding 
the future divergency problem in the CPL parametrization, and use them to probe the dynamics of dark energy 
in the whole evolutional^ history. Our divergency-free parametrizations are proven to be very successful in 
exploring the dynamical evolution of dark energy and have powerful prediction capability for the ultimate fate 
of the universe. Constraining the CPL model and the new models with the current observational data, we show 
that the new models are more favored. The features and the predictions for the future evolution in the new 
models are discussed in detail. 

PACS numbers: 95.36.+X, 98.80.Es, 98.80.-k 



The properties of dark energy are mainly characterized by 
the equation of state parameter (EOS), w. Extracting the in- 
formation of EOS of dark energy from observational data is 
very challenging owing to the accuracy of current data and 
our ignorance of dark energy. For probing the dynamical evo- 
lution of dark energy, under such circumstance, one has to 
parameterize w empirically, usually using two or more free 
parameters. Among all the parametrization forms of EOS, 
the Chevallier-Polarski-Linder (CPL) model [1] is the most 
widely used one and has been explored extensively. The form 
of the CPL parametrization is 

z 

w(z)-wo+wi- , (1) 

1 +z 

where z is the redshift, wq is the present-day value of the 
EOS, and W[ is the derivative of the EOS with respect to 
the scale factor a. The direct motivation of proposing such 
a parametrization form is to overcome the divergency feature 
of the linear form w{z) - wq + wizal high redshifts. Further- 
more, as Linder [1] suggested, the CPL parametrization has 
several advantages, such as a manageable two-dimensional 
phase space, well behaved and bounded behavior for high red- 
shifts, high accuracy in reconstructing many scalar field equa- 
tions of state, simple physical interpretation, etc. 

However, we have to point out that there exists a problem in 
the CPL model. The CPL model only explores the past expan- 
sion history properly, but cannot describe the future evolution 
due to the fact that \w{z)\ grows increasingly and finally en- 
counters divergency as z approaches - 1 . Undoubtedly, this is 
a nonphysical feature. Such a divergency problem prevents 
the CPL parametrization from genuinely covering the scalar- 
field models as well as other theoretical models. 

To overcome the shortcoming of the CPL model, we are 
interested in proposing novel parametrization forms of w(z). 
The new parametrizations will be contrived to exceed the CPL 
model entirely: inheriting the advantages of the CPL model, 
avoiding the disastrous divergency in the far future, and being 
more favored by the observational data. In this Letter, we are 
devoted to exploring more insightful parametrization forms 
for dark energy and probing the dynamics of dark energy in 
light of the novel parametrizations. 



The leading proposal we put forth for the EOS of dark en- 
ergy is of the form: 



w{z) - Wo + wi\ — ; In 2 

l+z 



(2) 



where wo also denotes the present-day value of w{z), and wi is 
another parameter characterizing the evolution of w{z). Note 
that a minus In 2 in the bracket is kept for maintaining wq to 
be the current value of w(z) and for easy comparison with the 
CPL model. Obviously, this new parametrization has well be- 
haved, bounded behavior for both high redshifts and negative 
redshifts. Thanks to the logarithm form in the parametriza- 
tion, a finite value for w(z) can be ensured, via the application 
of the L'Hospital's rule, in both limiting cases, z ^ oo and 
z ^ -1. This is the reason why we introduce a logarithm 
form in the new parametrization. For clearness, we list the 
values of w(z) in the limiting cases: 



w(z) 



Wo, 



Wo — wi In 2, 



forz = 0, 



forz 



-Hoo 



(3) 



Wo + wi(l - ln2), forz ^ -1. 



At low redshifts, the new form reduces to the linear one, 
w(z) ~ wq + wiz, where wi - -(ln2)wi. Of course, one 
can also recast it at low redshifts as the CPL form, w(z) ~ 
Wo + wiz/(l + z), where wi = (1/2 - ln2)wi. Therefore, 
it is clear to see that the new parametrization exhibits well- 
behaved feature for the dynamical evolution of dark energy. 
Without doubt, such a two-parameter form of EOS can gen- 
uinely cover many scalar-field models (including quintom 
models with two scalar fields and/or with one field with high 
derivatives) as well as other theoretical scenarios. 

It is well-known that dark energy drives the cosmic acceler- 
ation only at the late times (z ~ 0.5), whereas at the early times 
dark energy can be totally neglected due to the extremely low 
density compared to the matter or radiation component. Thus, 
one can only justify that the EOS of dark energy is around 
-1 in the recent epoch, but can tolerate more possibilities for 



2 



the early-time EOS of dark energy. For example, the EOS of 
dark energy might exhibit oscillating feature during the evolu- 
tion [2]. This is a fascinating possibility, deserving a detailed 
investigation. Based on this consideration, we further extend 
the above new parametrization (2) to an oscillating one, by 
replacing the logarithm function with a sine function: 



/ sin(l + z) 
w{z) = Wo -H wi I — — sin(l)). 



(4) 



Such a replacement is rather reasonable, lying in the fact that 
the two parametrizations roughly coincide in the recent epoch 
(low redshifts), since sin(l) w In 2 and cos(l) ^1/2 [note that 
sin(l) * 0.841, ln2 ^ 0.693, and cos(l) * 0.540]. Hence, the 
parametrization (4) describes the same behavior as the loga- 
rithm form (2) at low redshifts, but exhibits oscillating feature 
from a long term point of view. Also, we list the values of 
w(z) in the following limiting cases: 



w(z) = 



Wo, 



Wo - wi sin(l). 



forz = 0, 



forz 



Wo -H wi(l - sin(l)), forz 



+ 00 



-1. 



(5) 



We find that the two parametrizations, (2) and (4), also 
roughly coincide in the limiting cases, z — > oo and z — > -1. 

In what follows, we shall explore the dynamical evolu- 
tion of dark energy via the CPL parametrization and the 
new parametrizations. For convenience, we call the new 
parametrizations the logarithm form and oscillating form, re- 
spectively, hereafter Since our aim is to probe the dynamics 
of dark energy, we should try to avoid other indirect factors 
weakening the observational limits on the EOS; thus we as- 
sume a flat universe, = 0, consistent with the inflationary 
cosmology. From the Friedmann equation, the Hubble expan- 
sion rate can be written as 

H(z) = Ho [q„,(1 + zf + Q.(l + z)^ + (1 - n,„ - n,)f(z)f^ , 

(6) 

where = Qy(lH-0.2271A^eff), withQ-, = 2.469x10"^/!"^ for 
7"cmb - 2.725 K, A^eff the effective number of neutrino species 
(in this Letter we take its standard value, 3.04 [3]), and f(z) - 
oxp[3£dz'(l+w(z'))Kl+z')]. 

For constraining w(z), we use the current observational data 
from the type la supemovae (SN), the baryon acoustic oscilla- 
tions (BAO), and the cosmic microwave background (CMB). 
Such a combination of data sets is the most widely used one, 
sufficiently satisfying our aim of testing the new parametriza- 
tions and making a comparison. Of course, one can also add 
other data sets such as gamma-ray bursts, H(z), and so on, but 
we feel that this is not necessary for our present aim and leave 
a more sophisticated analysis to a future work with different 
goal. 

For the SN data, we use the 557 Union2 data compiled in 
Ref. [4]. The theoretical distance modulus is defined as 



Hthizi) = 5 log[(, DiiZi) + iio. 



(7) 



where jUo = 42.38 - 5 logjy h with h the Hubble constant Hq 
in units of 100 km/s/Mpc, and the Hubble-free luminosity dis- 
tance 



Dl{z)^{\+z) 



dz' 



E(z';ey 



(8) 



where E = H/Hq, and 6 denotes the model parameters. Cor- 
respondingly, the function for the 557 Union2 SN data is 
given by 



iMohs(Zi) - fi,h(Zi)f 

(rHzd 



(9) 



where cr is the corresponding Icr error of distance modulus 
for each supernova. The parameter fio is a nuisance parameter 
but it is independent of the data points. Following Ref. [5], 
the minimization with respect to fiQ can be made trivially by 
expanding of Eq. (7) with respect to fio as 



xlr,(0)^A-2^oB + filC, 



(10) 



where 



557 



A(0) = ^ 



/=i 



iMohsizi) - l^thizi; Ho = 0,d)f 
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HobsiZi) - firhizr, jUO = 0, 0) 



557 



1 



<Jzd 



i=l fob! 

Evidently, Eq. (10) has a minimum for /uo = B/C at 

B{0f 



xIj,{0)^A{0)- 



C 



(11) 



Since XsN,mw = ^SiV.mm' instead minimizing we will min- 
imize which is independent of the nuisance parameter yuo- 
For the BAO measurement, we use the data from SDSS 
DR7 [6]. The distance ratio (d-) at z = 0.2 and z = 0.35 
are 



^0.2 



rs(Zd) 
Dv{Q.2y 



do 



35 



rs(zd) 
Dy(0.35)' 



(12) 



where rs(zd) is the comoving sound horizon at the baryon drag 
epoch [7], and 



Dviz) = 



(X H(z')) 



Hiz) 



1/3 



(13) 



encodes the visual distortion of a spherical object due to the 
non Euclidianity of a FRW spacetime. The inverse covariance 
matrix of BAO is 



30124 -17227 
-17227 86977 



(14) 
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The function of the BAO data is constructed as: 



2 

Xbao 



(15) 



where c/, - (do2, do.a) is a vector, and the BAO data we use 
are c/o.2 = 0.1905 and do35 = 0.1097. 

The CMB is sensitive to the distance to the decoupling 
epoch via the locations of peaks and troughs of the acoustic 
oscillations. In this Letter, we employ the "WMAP distance 
priors" given by the seven-year WMAP observations [3]. This 
includes the "acoustic scale" U, the "shift parameter" R, and 
the redshift of the decoupling epoch of photons z,. The acous- 
tic scale Ia describes the distance ratio 0^(1,)/ r^iz*), defined 
as 



Ia = (1 



rsiz,) ' 



where a factor of (1h-z«) arises because Da(z*) is the proper an- 
gular diameter distance, whereas rsiz*) is the comoving sound 
horizon at z*. The fitting formula of r^iz) is given by 



rsiz) 



Jo 



da 



a^H(a) ^1 + (3Q/,/4f2y)a 



(17) 



where Q.b and Qy are the present-day baryon and photon den- 
sity parameters, respectively. In this Letter, we fix Qy = 
2.469 X 10"^/?"- and Q/, = 0.02246/;"- given by the seven- 
year WMAP observations [3]. We use the fitting function of 
z, proposed by Hu and Sugiyama [8]: 

z, = 1048[1 +0.00124(Qi/i2)-o"8][l +^i(n,„/i2)«2]^ (18) 



where 



2\ -0.238 



0.0783(Qfc/i^) 

1 -H39.5(Qi/l2)0.763 ' 



0.560 



1 -(-21.1(Qi/l2)l 



The shift parameter R is responsible for the distance ratio 
Da{z,)/H-\z,), given by [9] 



R(z.)= ^n,„Hl(l + z.)Da{z.). 



Following Ref. [3], we use the prescription for using the 
WMAP distance priors. Thus, the x~ function for the CMB 
data is 



XcMB 



(4 



where jc; - {Ia,R,z*) is a vector, and (C^ji^^)/; is the inverse 
CO variance matrix. The seven-year WMAP observations [3] 
give the maximum likelihood values: Ia{z*) - 302.09, 7?(zt) - 
1.725, and z* - 1091.3. The inverse covariance matrix is also 
given in Ref. [3]: 



2.305 29.698 
29.698 6825.27 
-1.333 -113.180 



-1.333 ^ 
-113.180 
3.414 



Since the SN, BAO and CMB are effectively independent 
measurements, we can combine the data sets by simply adding 
together the;if2 functions. Thus, we have 



2 _ ~2 , 2 ,2 

X - XsN Xbao + Xcmb ■ 



Note that^^^ is free of h, yN\a\e x bao ^^^Xcmb ^'■^^^ tqIq- 
vant to h. 

The three models have the same free model parameters, 
namely, - {0„,, wo, w\, h]. According to the joint data 
analysis, we obtain the best-fit parameters and the correspond- 
ing A',^,,,,- The best-fit, lcrand2crvaluesof the parameters with 
X^nin °f three models are all presented in Table 1. 

Figure 1 shows the likelihood contours for the three models 
in the wq - w\ and Q„, - h planes. For the CPL parametriza- 
tion, we have Q.,„ - 0.279, wq - -1.066, \v\ - 0.261 and 
h - 0.699, with y-. =544.186. We plot the likelihood con- 

' nun ^ 

tours for the CPL model in the panels (a) and (b) of Fig. 1. 
For the logarithm parametrization, the fitting results are Q.„, - 
0.280, wo = -1.067, wi = -1.049 and h = 0.697, with 
(16) xl,in = 544.081, smaller than that of the CPL model. The like- 
lihood contours for this case are shown in the panels (c) and 
(d) of Fig. 1. For the oscillating parametrization, we obtain 
the fitting results: Q,„ = 0.280, wo = -1.061, wi = -0.410 
and h - 0.695, with xl,i„ - 543.986, which is the smallest 
of the three. The likelihood contours for this case are shown 
in the panels (e) and (f) of Fig. 1. According to \ht x^„i„^ the 
new parametrizations are indicated to be more favored by the 
observational data. 

Next, we reconstruct the expansion history of the universe 
in light of the above fitting results for the three models, and 
then make a comparison for them. Since we focus on the prop- 
erties of dark energy, we only reconstruct the evolutionary be- 
haviors of the EOS of dark energy, w{z), and the deceleration 
parameter of the universe, q{z). The reconstructing results are 
shown in Fig. 2. As has been pointed out in the above, the 
CPL model has a problem: w{z) diverges when z approaches 
- 1 . So, the CPL parametrization can only properly describe 
the past evolution history but cannot genuinely depict the fu- 
ture evolution; it is incomplete in describing the evolutionary 
history of dark energy. Consequently, the CPL parametriza- 
tion is not capable of covering other dark-energy theoretical 
models. Such a problem can be explicitly seen in the recon- 
structed evolution plots, panels (a) and (b) of Fig. 2. We can 
see from these two plots, w{z) and q{z), that although the CPL 
model can do a good job in describing the past evolution of 
dark energy, it totally loses the prediction capability for the 
future evolution of dark energy. The novel parametrizations 
successfully overcome the shortcoming of the CPL model. 
The reconstructed evolutionary plots, panels (c) and (d) for 
the logarithm form (2) and panels (e) and (f) for the oscil- 
lating form (4), indicate that the both new models can nicely 
describe the whole evolution history of dark energy. 

Comparing the panels (a), (c) and (e) of Fig. 2, we find that 
all the three models favor a quintom behavior [10] that the 
EOS crosses -1 around the recent epoch. This w - -I cross- 
ing feature is only mildly favored by the CPL and the oscillat- 
ing models, but is explicitly favored by the logarithm model. 
For the CPL model, since it forfeits the prediction capability 
for the future evolution, we cannot say anything about the ul- 
timate fate of the universe. For the oscillating model, the fate 
of the universe is not definitely determined: the big rip may or 
may not occur For the logarithm model, the universe will def- 
initely move towards its tragic destiny: the cosmic doomsday 



(19) 



(20) 



(21) 



(22) 



(23) 
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TABLE I: The fitting values for the CPL, logarithm (Log) and oscillating (Sin) models. 



Model 




n„, 




Wo 








h 




A. mm 


CPL 


n 970+0.032 


.+0.047 
'•^'-'^•'-0.038 


(20-) 




(20-) 


0.261!^;^"^ (I'^)!i608 (2tr) 


699+" 

"•"^=^-0.034 


/.N +0.041 


(20-) 


544.186 


Log 




.+0.047 
l.^ "^^-0.040 


(2a-) 


-1.067-^(1^)-- 


(20-) 


-1.049!^-™^ (1^)-;™ (2^) 


697+0.031 


^^■^^-0.036 


(20-) 


544.081 


Sin 


0.280!»:«i 


r 1 ^\+0.050 
(l'^)-0.040 


(2cr) 


-1.061-1- (1^)!;;:^^^ 


(20-) 


-0.041!j-]^J(l^)!^-»(2<r) 


695+""^^ 

"•"^-'-0.038 


/ 1 „\+0.046 
(lo^)-0.054 


(20-) 


543.986 



CPL 
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FIG. 1: The probability contours at lo- and 2o- confidence levels in the wo - wi and £l„- 
(Sin) models. 



- h planes for the CPL, logarithm (Log) and oscillating 
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FIG. 2: The reconstructed evolutionary histories (from past to future) for w(z) and q(z) in the tliree models. 



will happen at about 3cr level. If we only focus on the past 
evolution, we find that the CPL model performs better than 
the oscillating model. The best one among the three is with- 
out doubt the logarithm model, not only for the description of 
the past evolution but also of the future evolution. From the 
panels (c) and (e) we see clearly that the oscillating model be- 
haves similarly to the logarithm within the range from z = - 1 
to z = 2, while the oscillating feature emerges in the oscillat- 
ing model from a long term perspective. Comparing the pan- 
els of (b), (d) and (f) of Fig. 2, we also find that the logarithm 



model performs the best. From the reconstructed q{z) plots 
we see that the accelerated expansion of the universe starts at 
a redshift around 0.5 - 0.7. According to the new models, for 
the future evolution, the change rate of the cosmic accelera- 
tion, \dq{z)/dz\, will first increase and then decrease, with the 
pivot around z ~ -0.45. We find that for the far future, the 
change rate \dq(z)/dz\ for the logarithm model is still rather 
large but for the oscillating model approaches zero. 

In summary, we have proposed two novel parametrizations 
for the EOS of dark energy, successfully avoiding the future 
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divergency problem of the CPL parametrization, and used 
them to probe the dynamics of dark energy not only in the 
past evolution but also in the future evolution. We pointed 
out that the CPL parametrization can only properly describe 
the past evolution history of dark energy but cannot genuinely 
depict the future evolution of dark energy owing to the diver- 
gency of w(z) as z approaches -1. Such a divergency fea- 
ture forces the CPL parametrization to lose its prediction ca- 
pability for the fate of the universe and to fail in providing a 
complete evolution history for the dark energy. Consequently, 
the CPL model cannot genuinely cover scalar-field models as 
well as other dark energy theoretical models. Our new pro- 
posals, the logarithm form w(z) - wq + wi(^^i^ - In 2) and 

the oscillating form w{z) - wq + wi( ^'"||T''"' - sin(l)), exhibit 
well-behaved features for the EOS of dark energy in all the 
evolution stages of the universe. We constrained the new mod- 
els by using the current observational data including the 557 
Union2 SN data, BAO data from SDSS DR7 and CMB data 
from 7-yr WMAP. The fitting results show that the oscillating 
parametrization gains a minimal among the three models, 
seemingly more favored by the data. However, by reconstruct- 



ing the whole evolutionary histories of w(z) and q(z) from past 
to future via the fitting results, we found that the logarithm 
parametrization is more tightly constrained by the data. The 
reconstruction results show that the w = - 1 crossing feature 
is favored. Furthermore, the new models predict that in the fu- 
ture the cosmic acceleration will first speed up and then slow 
down. In the logarithm model, the cosmic doomsday seems 
inevitable, which will happen at about 3cr level. We believe 
that the novel parametrizations proposed in the present Letter 
deserve further investigations. 
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